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We present results for the cross-correlation between the WMAP lst-year cosmic microwave 
background (CMB) anisotropy data and optical galaxy surveys: the APM and SDSS DR1 
catalogs. Our measurement of a positive CMB-galaxy correlation on large angles (9 > 4°) 
yields significant detections of the Integrated Sachs- Wolfe (ISW) effect and provides a new 
estimate of dark-energy in the universe, Qa ~ 0.69 — 0.86 (2a range). In addition, the 
correlated signal on small angles (9 < 1°) reveals the imprint left by hot intra-cluster gas in 
the CMB photons: the thermal Sunyaev-Zeldovich (SZ) effect. 



1 Introduction 



The pattern of primary temperature anisotropics of the cosmic microwave background (CMB) 
is expected to be distorted by the large-scale structures of the universe as microwave photons 
travel from the last scattering surface (z ~ 1100) to us (see e.gHS). On large-scales (i.e, angular 
scales comparable to large clusters and super-clusters as seen in projection) , the observationally 
favored flat ACDM models predict such distortion is mainly produced by the energy injection 
photons experience as they cross time-evo lving d ark-matter gravitational potential wells: the so- 
called integrated Sachs- Wolfe effect (ISW) ' ' I On smaller scales, the primary CMB anisotropy 
pattern is altered when primordial photons scatter off free electrons in the hot intra-cluster gas: 
the thermal Sunyaev-Zeldovich effect (SZ)^l. 

Although these secondary anisotropics can be potentially measured from CMB maps alone, 
in practice the ISW detection is severely limited by primary CMB anisotropies and cosmic 
variance, whereas the SZ requires high-spatial resolution, multi-frequency observations and it is 
potentially contaminated by point sources. Alternatively, large-scale structure tracers, such as 
galaxy surveys, provide a unique window to probe the baryon and dark-matter distribution at 



intermediate redshifts (z ;$ 2) and its imprint in microwave photons without significant confusion 
from other cosmological signal d ' • \ The combination of nearly full-sky high-sensitivity CMB 
maps obtained by WMAP^ with wide large-scale structure surveys has recently led to the first 
detec tions of the ISW effect, setting new constraints on the dark-ene rgy cont ent of the universe 
it h as a i so allowed probing the SZ effect with cluster templates ^^^ . 
Here we concentrate on the cross-correlation of the cosmic CMB anisotropics measured by 
WMAP with galaxy number count fluctuations in the APM Survey^ and the first data release 
of the Sloan Digital Sky Survey (SDSS DR1). Our CMB-galaxy correlation analyses find 
significant detections for both the integrated Sachs- Wolfe (ISW) and thermal Sunyaev-Zeldovich 
(SZ) effects. The reported ISW detection is in good agreement with previous analyses b ased 
on X-ray and radio sources^*^. Further details on t he r esults presented here are given in^^l. 
A similar analysis using SDSS data was presented in ^ while, more recently, ^ computed the 
CMB-galaxy correlation for 2MASS galaxies. 



2 Data 



We make use of the largest datasets currently available to study the CMB-galaxy cross-correlation. 
In order to probe the galaxy distribution, we have used the APM survey as well as selected sub- 
samples from the SDSS DR1. 

The APM Galaxy Survey is based on 185 UK IIIA-J Schmidt photographic plates each 
corresponding to 5.8 x 5.8 deg 2 on the sky limited to bj ~ 20.5 and having a mean depth of 
~ 400 Mpc/h for b < —40 deg and 5 < —20 deg. These fields where scanned by the APM 
machine and carefully matched using the 5.8 x 0.8 deg 2 plate overlaps. Out of the APM Survey 
we considered a 17 < bj < 20 magnitude slice, which includes 1.2 million galaxies at a mean 
redshift z = 0.15, in an equal-area projection pixel map with a resolution of 3.5', that covers 
over 4300 deg 2 around the SGC. 

The SDSS DR1 survey covers ~ 2000 deg 2 (i.e, 5 % of the sky). The samples we analyzed 
have different redshift distributions and a large number of galaxies (10 5 -10 6 , depending on the 
sample). We concentrate our analysis on the North sky (~ 1500 deg 2 , ie, 3.6 % of the sky), 
because it contains the largest and wider strips. The South SDSS DR1 (~ 500 deg 2 ) consists 
of 3 narrow and disjoint 2.5° strips, which are less adequate for our analysis. Our main sample, 
hereafter SDSS all, includes all objects classified as galaxies with extinction corrected magnitude 
r < 21, and a low associated error (< 20%). This sample contains ~ 5 million galaxies distributed 
over the North sky. Its predicted redshift distribution is broad and has a median redshift z ~ 0.3. 
Our color selected high-redshift sample (SDSS high-z thereafter) comprises ~ 3 x 10 5 galaxies, 
with z ~ 0.5. 

For the CMB data, we use the first-year full-sky WMAP maps 1 ^. We shall focus on the 
V-band (~ 61 GHz) as it conveniently combines low pixel noise and high spatial resolution, 21'. 
In addition, we have also used the W-band and a foreground "cleaned" WMAP mapl^to check 
that our results are free from galactic contamination. We mask out pixels using the conservative 
KpO mask, that cuts out 21.4% of the skyEI. All the maps used have been digitized into 7' pixels 
using HEALPix 



3 Cross-Correlation Estimation 



We define the cross-correlation function as the expectation value of density fluctuations 5q = 
Nq/ < Ng > — 1 and temperature anisotropies At = T — Tq (in fiK) at two positions fii and 
n2 in the sky: wtg(&) = {^-T (^1)^0(^2)}, where 9 = |na — fii|. 



"http:/ /www.sdss.org/drl 

6 http:/ /www.eso.org/science/healpix 
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Figure 1: (Left) error estimation for the WMAP-APM; (Right) errors for the WMAP-SDSS. Solid line gives 
errors from the dispersion in 200 Monte-Carlo simulations; Squares with errobars give the mean and dispersion 
in the jack-knife error estimation over the same simulations. Dashed lines show the jack-knife errors in the 
WMAP-Galaxy correlation for the real samples. Angular scales are given in degrees. 



We compute the CMB-galaxy correlation and the associated statistical error-bars using the 
jack-knife (JK) method (see^ and references therein for further details) The survey is divided 
into M separate regions the sky, each of equal area (M = 8 — 16 depending on the survey used; 
we find no significant differences as a function of M). The wtg analysis is then performed M 
times, each time removing a different region, the so-called JK subsamples. The covariance Cjj 
for wtg between scales 9i and 9j is obtained by re-scaling the covariance of the JK subsamples by 
a factor M — 1 To test the JK errors and their covariance we have also run 200 WMAP V-band 
Monte-Carlo (MC) realizations of the measured WMAP temper atu re angular power-spectrum^ 
and the white noise for the V-band estimated for 1 year of dataJEl. 

For each MC simulation we estimate the mean "accidental" correlation wtg oi simulated 
CMB maps to the SDSS galaxy density fluctuation map. We also estimate the associated JK 
error in each MC simulation. Fig. ^ compares error estimates for the APM (left panel) and 
SDSS (right) surveys. They show the 'true' sampling error from the dispersion of wtg{@) in 200 
MC simulations, along with the mean and dispersion of the JK errors in the same simulations. 
The JK and MC errors agree within 10 — 20% for 9 £ 5 degrees, i.e, the scales relevant for the 
ISW effect (see below). The observed differences in error estimation, particularly on smaller 
scales, are not totally surprising as the MC simulations do not include any physical correlations 
and use a CMB power spectrum that is valid for the whole sky, and not constraint as to match 
the CMB power over the region covered by the galaxy survey. Actually, the JK errors provide 
a model free estimation that is only subject to moderate (20%) uncertainty, while MC errors 
depend crucially on the model assumptions that go into the simulations. 

FigEl shows the CMB-galaxy correlation wtg for the APM (left) and SDSS (right) catalogs 
together with the corresponding JK error. We derive the significance of the detected correlation 
taking into account the large (JK) covariance between neighboring (logarithmic) angular bins 
in survey sub-samples. To assign a conservative significance for the detection (i.e. against 
wtg = 0) we estimate the minimum x 2 fit f° r a constant wtg an d give the difference A% 2 to 
the wtg = null detection. For example, at scales 9 = 4 — 10° we get: wtg = 0.35±0.13/iK for 
the APM, wtg = 0.26 ± 0.13^K for the SDSS all and w T g = 0.53 ± 0.21/iK for the SDSS high-z 
sample, in all cases we give 1-a errorbars. We find the largest significance in the CMB-galaxy 
correlation for the SDSS high-z sample: A% 2 = 9.1 (ie probability, P = 0.3% of no detection) for 
9 < 10°. In order to derive the significance levels for the ISW and SZ effects from the observed 
CMB-Galaxy correlations, we shall first discuss model predictions. 




Figure 2: (Left) WMAP-APM correlation: the two solid lines and the dotted line show wtg results for WMAP 
bands V, W, and the foreground "cleaned" map. Boxes show the 68% confidence levels. Also shown are the 
theoretical predictions for ISW and SZ (upper and lower short-dashed lines) and their sum (long-dashed line) 
for the best-fit ACDM model. (Right) WMAP-SDSS correlation: long dashed-line shows the measurement for 
the SDSS HIGH-Z sample, while the solid line displays the correlation for the SDSS ALL sample. For reference, 
short-dashed line displays the same measurement using the APM galaxy survey (see left panel) instead of SDSS. 



4 Comparison to Theoretical Predictions 

Photons coming to the observer from a given direction in the sky n suffer an ISW temperature 
change given by: AT /sw (n) = —2 J dz<&(h,z), and for a flat universe V 2 $ = — AirGa 2 p m 5. In 
Fourier space it reads, z) = —3/2Q m (Ho/k) 2 5(k,z)/a, and thus: 

vffifV) =< A™ W S G >=Jf P(k) g(kO) (1) 

being, g{k6) = 1/2-n J dz Wisw(z) Wq{z) jo(kOr), where the ISW window function is given by 
Wisw = -3£l m {H /c) 2 F(z), with c/H ~ 3000fr Mpc -1 , F = d{D/a)/dr = (H/c)D{f-l), and 
/ ~ quantifies the time evolution of the gravitational potential. The galaxy window 

function is Wq — b(z) D(z) 4>g(z), which depends on the galaxy bias, linear dark- matter growth 
and the galaxy selection function. The ISW predictions for the 3 samples are shown in in bottom 
panel of FigEl Unless stated otherwise, we use the concordance ACDM model with Q m = 0.3, 
tt A = 0.7, T ~ hn m = 0.2 and a 8 = 1. 

For the thermal Sunyaev-Zeldovich (SZ) effect, we assume that the gas pressure 5 gas fluctu- 
ations are traced by the galaxy fluctuations S gas ~ b gas 5g with a relative amplitude given by 
the gas bias, b gas ~ 2, representative of galaxy clusters. A simple estimate of the SZ effect is 

thusEDI, 

w sz {9) = _ bgas AT Wgg{0) (2) 

where AT = j(x)y To, is the mean temperature change in CMB photons Compton scattered 
by electrons in hot intracluster gas, To ~ 2.73K is the mean CMB temperature, y ~ 1 x 10~ 6 
is the mean Compton parameter induced by galaxy clusters estimated for our samples, and 
j(x) = —4.94 is the negative SZ spectral factor for the V-band. The weak lensing effect prediction 
is quite similar to the ISW, we just need to replace the time derivative of the Newtonian potential 
by its 2D Laplacian^l W Lens = 3k 2 n m (H /c) 2 (D/a)/d{r), d{r) being the angular distance to 




Figure 3: (Left) Theoretical predictions: (Bottom panel) Continuous, long and short dashed lines show the 
ISW, SZ and lensing predictions. Different sets of lines correspond to the APM (black), SDSS all (red) and SDSS 
high-z (blue) samples. (Top panel) The total prediction (ISW+SZ+Lensing) for the 3 samples. We have assumed 
a ACDM model with a fixed b gas = 2 in all cases, b = 3 for SDSS high-z and b = 1 for APM and SDSS all. 
(Right) Predictions and measurements as a function of redshift: squares and circles correspond to measured 
cross-correlation at 9 — 0.6° and 6 — 6° respectively. Continuous lines show the corresponding SZ+ISW+lens 
predictions for b = 1 and b gas = 2. Dotted and dashed lines shows the ISW with b = 1 and b — 3 . 



the lensing sources. The total predicted correlation is thus the sum of three terms: the ISW, 
thermal SZ and Lensing contributions, wtg = w t S <3^ + w tg + w tg & • FigEl shows individual 
contributions of these effects (bottom panel) and the total (top) for the 3 samples analyzed. The 
ISW effect typically dominates for angles 9 > 4°, while the SZ effect is expected to be significant 
on small scales 9 < 1°. Lensing is found to be negligible at all scales for our samples. 

Before one can make a direct comparison between theory and observations, the issue of 
galaxy bias must be addressed. As shown in Fig |3] (see plot on the right), the higher (SDSS) 
redshift sample requires a high bias b > 1 to explain the large cross-correlation seen at all scales. 
At low redshifts the measured correlation is dominated by the thermal SZ on small scales 9 < 1°, 
while the ISW still dominates on large-scales, 9 > 4°. Here no bias is required to reproduce 
the observations. Actually, this agrees quite well with our self-consistent bias estimation: for 
each sample we estimate the ratio b 2 = wgg I W M M , where wmm and wgg are the (theoretically 
predicted) matter and (measured) galaxy auto-correlation functions. For APM and SDSS all 
samples we find b 2 ~ 1 , while for the SDSS high-z sample we get b 2 ~ 6. 

4-1 Significance Tests 

ISW effect: On large scales 9 > 4°, the ISW effect is expected to dominate for all survey 
depths (see Fig EJ) • Therefore, from the large-angle CMB-Galaxy correlation, we can directly 
infer the ISW effect (ie, wtg = w t<3 > see enc ^ °f ©• ^ n particular, for the APM survey, a 
constant correlation fit rejects the null detection with high significance Ax 2 = 6.0 (P = 1.4%), 
similar to the values for the SDSS high-z sample Ax 2 = 6.1 (P = 1.3%). A smaller significance 
is obtained from the SDSS all sources: Ax 2 = 3.9 (P= 4.8%). Since these samples are basically 
independent, we can combine them to infer a total significance for the ISW detection: we find 
a total Ax 2 = 16 (P = 0.1% for 3 d.o.f) corresponding to a 3.3o\ Note we could do better 
using a (scale dependent) ACDM model theory prediction, but at the cost of introducing model 
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Figure 4: Estimating dark-energy: Long-dashed, short-dashed and dot-dashed lines show the probability distri- 
bution for £2a in the SDSS all, APM and SDSS high-z samples. Solid line shows the combined distribution. 

dependent detection levels. Moreover, we can further include the ISW-dominated small angle 
bins in our deepest sample, where SZ is negligible, increasing the significance to A% 2 = 18.8, 
(P = 0.03% for 3 d.of.), ie we detect the ISW effect at the a 3.6<r level. 

SZ effect: We can estimate the significance of the drop in the signal at small angles in the 
APM and SDSS all samples due to the SZ effect (see Fi^2l using the best-fit constant at large 
angles (ie, the ISW signal) and ask for the observed deviation from such value at smaller scales. 
For 6 < 1°, we find: tofg = -0.41 ± 0.16 for APM, and iof§ = -0.27 ±0.11 for SDSS all (1-ct 
errorbars). This test gives Ax 2 = 8.5 (P = 0.3%) for the APM sample and A X 2 = 5.5 (P = 2%) 
for the SDSS all. 




5 Discussion 

We have cross-correlated WMAP with the APM and SDSS DR1 optical galaxy surveys. Our 
analysis includes 5800 deg 2 (14 % of the sky) and comprises 7.6 million galaxies in compact 
regions of the north and south hemispheres. We obtain significant cross-correlations for galaxy 
samples in a wide redshift range (z ~ 0.15 — 0.5). We detect a positive large-scale correlation 
which is in good agreement with the ISW effect and its redshift evolution as expected from 
ACDM models. The combined analysis for the 3 samples investigated (APM, SDSS-all and 
SDSS high-z) yields a 99.97% ISW detection level (3.6a). 

Before using matter predictions to estimate parameters from the observed CMB-galaxy 
correlation, we self-consistently estimated the galaxy bias b by comparing the matter angular 
auto-correlation function in each model to the measured galaxy auto-correlation in each sample. 

Once galaxy bias is determined, we find that all the measured cross-correlations on large 
scales are in good agreement with ISW predictions for a dark-energy dominated universe. Fig0] 
shows the probability distribution for £l\ in a flat ACDM model. We have fixed as = 1, h = 0.7 
and VLm + = 1- As we vary £l\ the shape parameter for the linear power spectrum P(k) 
consistently changes T = MIm ■ We only use the data for 8 > 4° , where the ISW is the dominant 
contribution for all samples. 

From our analysis a coherent dark-energy dominated universe arises: both APM and SDSS 
samples point to large values of with the best fit £l\ ~ 0.8 and a rather narrow 2a range 

= 0.69 — 0.87. We stress that this dark-energy estimation is independent from other known 
probes (e.g, SN type la data). 



We have also found evidence for the thermal SZ effect from the drop of the CMB-galaxy 
correlation on small-scales in the low-z samples of APM and SDSS galaxies. The estimated SZ 
effect is compatible with a Compton parameter y ~ 1 x 10~ 6 . These new measurements can be 
used to constrain the redshift evolution of the physical properties of gas inside galaxy clusters. 

Upcoming wider and deeper Galaxy surveys (e.g future data releases from the SDSS) in 
combination with higher-spatial resolution and sensitivity full-sky CMB maps (e.g. 2yr-WMAP 
data, PLANCK) shall eventually allow for a higher significance detection of both the ISW 
and SZ signals. Moreover th ese new g eneration surveys can potentially extract the challenging 
lensing signal on large-scales I 23 I 12 U4I This program will eventually provide us with a better 
understanding of the gravitational instability picture and set new and tighter constraints on 
basic cosmological parameters. 
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